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ABSTRACT 

Linear  atability  theory  ia  being  used  to  inveatigate  the  mechanisms 
which  lead  to  boundary  layer  transition  on  re-entry  vehicle  noBetipa. 

Specific  efforts  are  aimed  at  determining  the  effects  of  various  parameters 
including  surface  roughness,  ablation,  wall  cooling,  local  Mach  numbers, 
pressure  gradient,  and  axisymmatric  geometry  on  transition,  toughness  is 
simulated  by  means  of  the  "turbulent  sublayer"  model  which  has  been  extended 
to  compressible  flow.  Results  indicate  that  both  roughness  and  ablation 
strongly  reduce  the  transition  Reynolds  number.  A series  of  parametric 
stability  solutions  has  been  calculated  for  non-zero  pressure  gradient  cases. 
Raaults  show  that  even  at  the  subsonic  Mach  numbers  which  characterize  nose- 
tip  boundary  layers,  the  effects  of  compressibility  have  substantially  altered 
the  stability  characteristics  of  the  boundary  layer.  Computations  of  the 
stability  characteristics  of  an  actual  nosetlp  boundary  layer  indicate  that 
for  smooth-wall  cases,  all  types  of  disturbances  are  stable  and  that  the 
margin  of  stability  is  considerable.  Inclusion  of  realijtic  roughnesses 
rapidly  lowers  the  critical  Reynolds  number.  Thus  far,  rough-wall  solutions 
have  been  completed  only  near  the  stagnation  point  where  the  boundary  layer 
ia  still  stable  to  all  disturbances.  Additional  computations  further  around 
the  noaetip  are  in  progress. 
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NOMENCLATURE 

Local  amplitude  of  disturbance  In  boundary  layer. 

Initial  amplitude  of  boundary -layer  disturbances. 

Constant  in  wall  roughness  model,  see  Eq.  (18). 

Dlmanslonlasa  blowing  rate,  B'  ■ (e^v^/p^u^Jst. 

Density-viscosity  product,  see  Eq.  (17). 

Specific  heat  at  constant  pressure. 

Velocity  variable  in  transformed  momentum  equation; 

F*  * u/u  . Also,  non-dimensional  frequency,  F ■mv/u *. 

Non-dimensional  frequency  in  rotated  coordinate  system; 

I*  • F/coe2iJ». 

Non-dimensional  anthalpy  in  tranaformed  energy  equation; 

G - K/H#. 

Stagnation  or  total  anthalpy;  also,  shape  factor,  H « 6*/6. 

Constant  in  wall  roughness  modal,  see  Eq.  (18). 

Molocular  thermal  conductivity;  also  used  to  signify 
roughness  height. 

"Total"  thermal  conductivity;  sum  of  molecular  plus 
"turbulent"  conductivity  in  near  wall  region. 

Mach  number. 

Pressure. 

Prandtl  number,  uC  /k. 

P 

"Total"  or  effective  Prandtl  number,  aea  Eq.  (9). 

Square  root  of  Reynolds  number,  R * Ae^. 

Reynolds  number  in  rotated  coordinate  system,  % - Rcosij*. 
Locus  radius  of  axisymraetric  body. 

Reynolds  number. 

Roughness  Reynolds  number,  Re^  ■ u^k/v^. 
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NOMENCLATURE  (cant’d) 

Arc-length  distance  along  surface  measured  from 
stagnation  point, 

. : J i , 

Stanton  number. 

■?; 

Temperature . 

TJ*a.  . '-'Iv 

Velocity  components  In  cartesian  (or  exisynnetric) 
coordinate  eye tea.' 

Velocity  at  top  of  roughnaaa  elaaent. 

Cartesian  coordinate  system;  x le  etreeawise  direction; 
y may  also  refer  to  radial  coordinate  axiaymmatric  ayetea. 

Complex  Wave  number  of  fluctuating  quantity. 

2,  dMe 

Yalkner-Skatn  pressure  gradient  pereneter.  0 - - j -Tjf  in 
Lavy-Laae  coordinates.  e 

Cone tent  In  wall  roughnaaa  model,  see  Eq.  (18). 

Ratio  of  specific  heats. 

Boundary-layer  thickness. 

Representative  boundary-layer  thickness. 

Boundary-Layer  displacement  thickness. 

Eddy  viscosity  for  momentum  transport. 

Eddy  conductivity  for  transport  of  haat. 

Transformed  normal  (to  the  well)  coordinate  in  Levy-Lees 
system;  eee  Eq.  (12). 

Momentum  thickness  of  boundary  layer. 

Molecular  viscosity. 

"Total"  viscosity,  sum  of  molecular  plus  "turbulent" 
viscosities  in  near-wall  region. 

Kinematic  viscosity,  v ■ \if p. 

Transformed  streamwisc  coordinate  in  Levy-Lees  system; 
see  Eq.  (11). 
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Superscript! 


NOMENCLATURE  (cont’d) 


Density. 

Skin  friction  at  the  wall. 

Refers  to  a general  variable;  azimuthal  coordinate  in 
cylindrical  coordinate  system. 

Angle  of  skewness  of  three-dimensional  wave.  *|»  - 0 
corresponds  to  wave  whose  direction  of  propagation  is 
in  the  x-direction. 

Frequency  of  fluctuating  quantity. 


Signifies  time-averaged  part  of  unsteady  quantity. 

Signifies  fluctuating  part  of  unsteady  quantity; 
also  indicates  differentiation. 

Signifies  amplitude  of  a fluctuating  quantity. 

Exponent  relating  viscosity  to  temperature. 


Subscripts 


Critical  condition. 

Conditions  evaluated  at  the  edge  of  the  boundary  layer. 

Imaginary  part  of  a complex  number. 

Conditions  evaluated  at  the  top  of  roughness  elements; 
or  referring  to  roughness. 

Conditions  evaluated  at  the  wall. 

Conditions  evaluated  far  upstream  from  the  body. 
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1.  INTRODUCTION 

As  a vehicle  re-enters  the  atmosphere,  its  entire  surface  is  initially 
embedded  within  an  all-laminar  boundary  layer,  and  transition  to  turbulence 
occurs  in  the  wake  behind  the  body.  During  the  descent,  the  location  of 
transition  moves  progressively  closer  to  the  body,  reaches  the  rear  skirt 
of  the  vehicle,  and  continues  to  move  toward  the  nose  of  the  vehicle  in  a 
more  or  less  steady  fashion.  During  this  time,  certain  anomalies  in  the 
flight  path,  or  in  atmospheric  conditions,  may  cause  the  transition  location 
to  remain  stationary  for  a time,  or  even  to  move  backwards  momentarily;  how- 
ever, after  Borne  critical  condition  has  been  surpassed,  the  transition  loca- 
tion ceases  this  gradual  forward  motion  and  instantaneously  jumps  from  the 
conical  section  of  the  vehicle  onto  the  noaetip.  Experiments  and  flight  data 
alike  have  shewn  that  the  nosetip  location  is  a stable  location  and  that 
transition  la  not  likely  to  jump  again  to  the  aft-end  of  the  vehicle.  After 
the  transition  point  has  jumped  to  the  nosetip,  it  once  again  resumes  its 
forward  progression,  moving  toward  the  noaetip  stagnation  point. 

When  the  transition  region  is  on  the  nosetip,  a potentially  strong  inter- 
action exists  between  the  high  local  heat  transfer  rates  which  occur  in  the 
vicinity  of  the  transition  region  and  the  regression  rate  of  the  ablative 
noaetip  surface.  Of  course,  the  surface  of  the  entire  nosetip  regresses  as 
the  body  descends,  but  the  high  transitional  heating  causes  increased  local 
ablation  rates  which  can  lead  to  radical  changes  in  the  shape  of  the  nosetip, 
especially  if  the  transition  location  remains  at  a fixed  position  on  the  nosetip 
for  a time.  A more  complete  understanding  of  the  mechanisms  which  cause  transi- 
tion to  occur  on  the  nosetip  would  enable  the  designer  to  minimize  the  thickness 
of  the  ablative  shield  which  is  used  on  the  nosetip,  while  still  maintaining  an 
adequate  margin  of  safety,  and  to  predict  more  accurately  the  ballistic  tra- 
jectory of  the  re-entering  vehicle.  Such  knowledge  could  also  be  used  to  select 
the  flight  path  in  a manner  which  would  minimize  the  tendency  to  create  radical 
shape  deformations. 

The  mechanisms  which  control  the  position  of  the  transition  region  when 
it  occurs  on  the  nosetip  are  many.  Among  them  are  the  free-stream  Mach  number, 
the  local  Mach  number  distribution  at  the  edge  of  the  boundary  layer,  the  Reynolds 
number,  the  wall  temperature,  the  surface  roughness  and  the  surface  ablation  rate. 
Of  these  various  parameters,  the  bluck  of  experimental  transition  data  which  has 
been  amassed  in  the  past  few  years  under  the  PANT  Program  (Powara,  1973;  Anderson, 
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1974)  suggests  that  the  roughness  characteristics  of  the  nosetip  surface  most 
strongly  Influence  the  location  of  transition.  A second  significant  parameter 
is  the  wall-to-frei-stream  enthalpy  ratio  which  can  have  particularly  Important 
effects  in  the  presence  of  wall  roughness.  Although  these  experimental  results 
have  greatly  increased  our  knowledge  of  the  transition  behavior  on  the  nosetip, 
significant  gaps  still  exist  in  our  understanding.  Effects  of  changes  in  body 
geometry  (either  design  changes  or  changes  caused  by  ablation)  are  poorly 
understood,  as  are  effects  of  scaling  the  body  slse.  This  report  summarises 
the  current  results  and  status  of  an  on-going  analytical  research  program  which 
is  aimed  at  obtaining  improved  understanding  of  the  mechanisms  which  lead  to 
transition  on  the  nosetip,  and  at  identifying  the  effect  of,  and  the  relative 
Importance  of,  the  various  parameters  which  affect  transition.  Although  the 
effects  of  roughness  ar-2  ablation  are  dealt  with  in  the  present  report,  only 
their  direct  effect  on  transition  is  considered;  no  description  of  the  coupling 
between  the  boundary  layer  and  the  shape  of  the  surface,  or  its  roughness  char- 
aceteristics,  is  included.  Finally,  it  is  noted  that  even  though  nosetipe  with 
rough  aurfaceo  and  finite  rates  of  heat  transfer  are  of  primary  lnterast,  the 
study  of  smooth-walled  nosetipe  and  adiabatic  nosetips  can  also  be  helpful  in 
elucidating  transition  meehanisoa.  Some  discussion  of  these  effects  is  pre- 
sented later. 

2.  TRANSITION  MECHANISMS  AND  THEIR  THEORETICAL  DESCRIPTION 

Before  discussing  methods  for  describing  the  transition  phenomena,  it  is 
desirable  to  outline,  in  aa  much  detail  as  possible,  the  events  which  occur 
during  boundary  layer  transition.  The  most  careful  and  thorough  experiments 
which  have  been  conducted  to  examine  transition  from  a "microscopic"  view- 
point are  the  series  of  experiments  which  have  been  carried  out  at  the  National 
Bureau  of  Standards  (Schubauer  and  Skramstsd,  1948;  Schubauer  and  Klebanoff, 
1956;  Klebanoff,  Tidstrom,  and  Sargent,  1962).  Although  these  measurements 
were  taken  in  incompressible,  flat  plate  boundary  layers,  we  anticipate  that 
many  of  the  same  phenomena  which  were  observed  in  these  experiments  will  also 
take  place  in  a typical  nosetip  boundary  layer.  To  be  sure,  certain  differences 
may  also  be  expected;  these  are  mentioned  later. 
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2.1  Experimental  Observations 

If  we  observe  the  development  of  a boundary  layer  from  the  leading  edge 
of  a flat  plate  to  a distance  which  is  sufficiently  far  downstream  that  the 
boundary  layer  has  become  fully  turbulent,  we  can  distinguish  a number  of 
distinct  regions.  In  the  first  region,  the  boundary  layer  is  laminar,  and 
if  small  disturbances  are  introduced  into  the  boundary  layer,  they  will  decay 
in  amplitude  aa  they  are  swept  downstream.  After  the  boundary  layer  has 
reschsd  a particular  thickness  (i.e.j  when  the  local  Reynolds  number  becomes 
large  enough),  disturbances  within  a particular  narrow  frequency  band  will  be 
amplified  as  they  are  swept  downstream,  while  others  will  continue  to  be  damped. 
Comparisons  with  boundary  layer  stability  theory  Indicate  that  those  waves  which 
grow  (as  wall  as  those  which  decay)  era  closely  described  by  the  linear  theory. 
Consequently,  this  region  can  be  referred  to  as  a region  of  linear  growth. 

Following  the  region  of  linear  growth  is  a region  of  non-linear  growth. 

The  precise  location  st  which  the  non-linear  effects  become  important  depends 
upon  the  amplitude  of  the  disturbances  lnsido  the  boundary  layer.  This  has  been 
demonstrated  experimentally  by  Klebsnoff  et  al.  (1962)  by  Introducing  artificial 
disturbances  into  the  boundary  layer  upstream  of  the  non-linear  growth  region 
and  observing  their  effect  on  the  transition  location.  As  the  amplitude  of 
these  external  disturbances  was  increased,  the  length  of  the  linear  region  de- 
creased, non-linear  effects  btgan  sooner,  and  the  transition  location  moved 
toward  the  leading  edge  of  the  plate.  For  sufficiently  large  external  distur- 
bances, the  linear  growth  region  was  completely  eliminated.  The  manner  in 
which  the  non-linearities  enter  is  particularly  significant.  They  first  appear 
as  a modification  of  the  mean  flow  profiles;  the  build-up  of  harmonic  frequen- 
cies in  this  initial  non-linear  region  is  observed  to  be  small.  The  physical 
meaning  of  this  interaction  between  the  unsteady  disturbances  and  the  mean  flow 
is  obvious:  energy  is  extracted  from  the  mean  flow  to  feed  the  growing  dis- 

turbances. In  the  simplest  theoretical  description,  this  interaction  is 
described  by  the  non-zero  correlation  between  the  streamwise  and  the  vertical 
velocity  fluctuations,  <u'v*>. 

One  remaining  aspect  of  the  non-linear  growth  region  in  flat  plate, 
transitional,  boundary  layers  is  of  particular  interest.  First  of  all,  these 
non-linear  effects,  even  in  a two-dimensional  boundary  layer,  havo  a pronounced 
three-dimensional  nature.  The  interaction  between  the  waves  and  the  mean  flow 
result  in  periodic  spanwise  variations  in  the  mean  flow  profiles.  (The 
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axisymmetric  realization  of  these  planar  three-dimensional  effects  would  appear 
u to  be  a significant  unknown  quantity  in  the  nosetlp  problem.)  At  some  spanwlse 

locations,  the  profile  has  become  slightly  fuller  than  the  Blasius  profile,  while 
at  intermediate  locations,  a distinct  point  of  inflection  has  been  developed.  It 
la  in  this  inflectional  profile  that  the  onset  of  the  next  region  is  noted, 

{j  namely,  the  intermittent  appearance  of  turbulent  spots.  These  turbulent  spots 

are  generated  in  the  inflectional  portions  of  the  boundary  layer  profile  be- 
cause of  a secondary  instability  which  exhibits  a very  rapid,  almost  explosive, 
amplification  which  causes  the  high  frequency  content  of  the  disturbances  to 
U Increase  rapidly,  resulting  in  the  birth  of  a turbulent  spot.  It  is  interesting 

to  note  that  this  streamwlse  location  (where  turbulent  spots  begin)  is  one  of  the 
earliest  locations  which  is  used  in  experimental  investigations  as  a measure  of 
the  "beginning"  of  transition.  However,  as  noted,  a variety  of  fundamentally 
' > important  transition  phenomena  have  Already  taken  place  before  the  turbulent 

spot  appears. 

Following  the  appearance  of  turbulent  spots,  comes  the  final  region  in 
which  turbulence  spreads  to  encompass  the  entire  boundary  layrr  leading  to 
the  familiar  fully  turbulent  characteristics.  In  this  region  the  mean  velocity 
profiles  develop  their  very  full  turbulent  shape,  and  the  skin  friction  and 
heat  transfer  coefficients  at  the  wall  increase  to  their  fully  turbulent  values. 

In  summary,  a number  of  observations  can  be  made.  First  of  ail,  it  is  clear 
that  significant  "transitional"  events  occur  before  the  "beginning"  of  transition 
as  measured  in  typical  experiments.  Any  method  which  is  aimed  at  determining  the 
mechanisms  which  lead  to  transition  must  Include  these  phenomena.  Second,  the 
initial  mean  velocity  profile  deformations  which  occur  in  the  boundary  layer  are 
three-dimensional  in  nature,  but  it  is  the  inflectional  profiles  which  act  to 
generate  turbulence.  As  the  profile  deforms  from  its  original  Blasius  shape  to 
its  final  turbulent  shape,  it  first  deforms  to  an  inflectional  profile  (£.«.,  a 
less  full  profile)  and  then  reverts  to  the  full  turbulent  profile.  Thus,  the 
deformation  of  the  velocity  profile  is  not  a mono tonic  deformation  from  the 
Blasius  shape  to  the  turbulent  profile.  A third  observation  is  that  the  initial 
linear  amplification  region  frequently  consumes  the  largest  percentage  of  the 
distance  to  the  initial  turbulence  breakdown  location.  In  other  words,  the  non- 
linear amplification  mechanisms  are  quite  strong  and  take  place  very  rapidly 
(over  a short  distance);  implicit  use  of  this  observation  has  been  made  in  our 
linear-stability  based  analysis  which  is  described  below.  Finally,  it  is  noted 
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that  the  waves  which  are  observed  to  gro-j  in  the  pre-transitional  boundary 
layer  are  much  lower  in  frequency  than  are  typical  turbulence  phenomena. 

The  higher  frequencies  which  are  characteristic  of  turbulence  do  not  appear 
until  the  secondary  instability  has  begun. 

All  the  above  discussion  pertains  to  transition  on  a flat  plate,  and 
although  we  can  expect  to  see  some  similarities  between  the  flat  plate  and 
the  nosetip  boundary  layers,  we  also  anticipate  that  there  will  be  some  dif- 
ferences. The  most  striking  difference  is  the  pressure  gradient  in  the  two 
problems.  Whereas  the  above-described  flat  plate  data  was  obtained  under 
constant  free-stream  velocity  and  sero  pressure  gradient  conditions,  the 
nosetip  flow  field  is  characterized  by  the  strongly  favorable  pressure 
gradient  conditions  which  are  representative  of  the  expansive  flow  away 
from  an  axisymmetrlc  stagnation  point.  (If  the  pressure  gradient  on  the 
nosetip  is  expressed  in  terms  of  the  familiar  Falkner-Skan  parameter,  ft, 
the  noeatlp  has  a value  of  0 - 0.5,  while  the  flat  plate  has  a value  of  0 -0.) 
One  might  expect  the  mechanisms  which  control  transition  in  this  highly  favor- 
able pressure  gradient  to  be  somewhat  different  than  those  which  are  observed 
on  a flat  plate.  Our  analysis  bears  this  out.  Similar  differences  might  be 
expected  due  to  the  axisynmetric  geometry,  the  compressibility  effects  and 
the  wall  cooling  effects  which  occur  on  the  nosetip.  Ablation  and  wall  rough- 
ness are  also  expected  to  Influence  the  transition  mechanisms  on  a nosetip  as 
indicated  above. 

2.2  Methods  for  Predicting  Transition  Locations 

A number  of  methods  for  predicting  the  location  of  transition  have  been 
proposed.  These  methods  range  from  purely  empirical  ones  to  methods  which 
solve  the  complete  unsteady  Navier-Stokes  equations.  Empirical  methods  are 
generally  concerned  with  correlating  the  macroscopic  results  of  transition 
testing;  that  is  to  say,  empirical  methods  try  to  predict  the  magnitude  and 
the  direction  of  the  change  in  the  transition  Reynolds  number  which  is  caused 
by  a corresponding  change  in  some  parameter  (such  as  the  Mach  number,  the  free- 
stream  unit  Reynolds  number,  eto. ).  Quite  naturally,  most  correlation  tech- 
niques are  based  partially  on  theoretical  principles  (and  as  such  should  be 
classified  semi-empirical  theories) . 

Some  of  the  more  recent  correlation  techniques  which  are  used  for  the 
nosetip  transition  problem  are  the  one  which  has  been  developed  by  Anderson 
(1973)  and  the  one  being  developed  by  White  (1974).  These  techniques,  which 
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combine  current  theoretical  understanding  of  the  nosetip  problem  with  the 
available  experimental  data,  have  resulted  In  useful  engineering  correlations 
for  a limited  range  of  parameters. 

The  extension  of  phenomenological  turbulence-modelling  theories  to  low 
Reynolds  number  ranges  in  the  hope  of  describing  the  characteristics  of  the 
transitional  boundary  layer  has  also  become  Increasingly  popular.  This 
approach  was  originally  proposed  by  Donaldson  (1969),  and  applications  of 
his  basic  approach  to  the  nosetip  problem  are  currently  being  attempted  by 
Vinson  st  at.  (1974)  and  Wilcox  (1974).  The  turbulence  modelling  approach, 
whan  combined  with  a finite-difference  solution  of  the  equations  of  motion, 
has  the  attractive  appeal  that  it  gives  "exact"  results  in  the  laminar  region 
of  the  boundary  layer,  and  contains  a phenomenological  turbulence  modal  which 
has  been  adjusted  to  give  reasonable  results  in  the  turbulent  flow  regime. 

Thus,  the  only  portion  of  the  boundary  layer  which  is  missing  is  the  transi- 
tional part.  Several  turbulenca  codes  (fi.g.,  Mellor  and  Herring,  1972;  Beckwith 
and  Bushnall,  1967)  include  this  region  by  a "switch"  which  gradually  "turns 
on"  the  turbulence  model  during  the  transitional  region.  The  location  of 
transition  is  specified  by  means  of  a semi-empirical  correlation.  No  attempt 
is  made  to  Include  any  of  the  physics  of  transition  in  the  computation  of  the 
transitional  region.  The  approaches  of  Donaldson,  Vinson  s t at. , and  Wilcox 
represents  attempts  to  extend  the  applicability  of  tha  turbulence  models  to 
low  Reynolds  numbers  by  including  some  of'  the  physics  of  transition  in  the 
model.  These  approaches,  then,  basically  become  experimental  data  correlation 
techniques  which  use  partial  differential  equations  rather  than  algebraic  equa- 
tions. Whether  or  not  any  attempt  is  made  to  model  some  of  the  physics  in  the 
transitional  boundary  layer,  turbulence-model  approaches  have  the  advantage 
(and  the  disadvantage])  that  they  always  yield  a solution  to  the  complete 
laminar-transitional-turbulent  boundary  layer.  Further,  if  sufficient  data 
exists  to  define  the  arbitrary  constants  and  functions  in  the  models,  these 
turbulence  approaches  can  be  effectively  used  to  Interpolate  between  existing 
sets  of  experimental  results.  Since  they  contain  only  a minimal  amount  of 
physicB  (particularly  as  regards  the  early  portions  of  transition  where  the 
turbulence  which  they  model  does  not  yet  exist),  the  turbulence  model  approaches 
cannot  be  reliably  extrapolated  to  flow  regimes  for  which  they  have  not  been  ad- 
justed, nor  can  they  be  used  to  study  the  mechanisms  of  transition.  Like  the 
other  data-correlation  methods,  they  can  only  attempt  to  predict  macroscopic 
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effects  in  the  transitional  boundary  layer.  A further  discussion  of  the  use 
of  turbulence  models  for  transition  calculations  and  comparison  with  the  linear 
stability  technique  is  given  in  Section  5. 

The  moat  ambitious  approach  to  predicting  boundary-layer  transition  is 
through  the  numerical  solution  of  the  full,  unsteady  Navier-StokeB  equations. 

It  la  normally  assumed  that  the  solutions  of  the  unsteady  Navier-Stokes 
aquations  contain  the  complete  phenomena  of  transition;  however,  as  indicated 
above,  the  full  threa-dlmsnalonal  equations  must  be  used  to  obtain  realistic 
results.  This  requires  enormous  amounts  of  computational  storage  and  time, 
and  problems  of  sub-grid  scale  fluctuations  must  still  be  overcome.  Despite 
the  difficulties,  progress  in  this  direction  has  been  reported  by  Orssag  (1974) 
and  Grosch  (1974),  but  at  the  present  time  this  technique  remains  too  expensive 
for  use  In  predicting  the  flows  over  realistic  bodies. 

2 . 3 The  Use  of  Linear  Stability  Theory  for  Investigating  the  Mechanisms 

of  Transition 

The  method  we  have  used  for  our  analysis  of  the  mechanisms  which  lead  to 
transition  is  based  on  linear  stability  theory.  Like  the  numerical  solution 
of  the  complete  Navier-Stokea  equations,  the  linear  stability  approach  repre- 
sents an  exact  solution  of  the  unsteady  Navler-Stokas  equations  in  the  early 
portions  of  the  pre-transltlonal  boundary  layer.  However,  like  the  seml- 
emplrlcal  correlation  techniques  and  the  turbulence  models,  the  linear  sta- 
bility approach  cannot  predict  transition  locations  without  soma  empirical 
input.  The  reason  for  this  is  simply  that  transition  is  wn  Inherently  non- 
linear phenomena.  For  example,  transition  cannot  occur  until  non-linear 
Interactions  between  the  wave-like  disturbances  and  the  mean  flow  profiles 
have  taken  place.  Nevertheless,  as  indicated  above,  much  of  the  controlling 
growth  of  disturbances  inside  the  boundary  layer  can  be  accurately  described 
by  the  linear  theory  in  many  instances,  and  this  represents  virtually  the 
only  promising  theoretical  approach  to  take  if  an  understanding  of  the  transi- 
tion mechanism  is  desired. 

The  equations  of  linear  stability  theory  are  obtained  from  the  unsteady, 
Navier-Stokes  equations  by  expressing  each  of  the  dependent  variables  in  terms 
of  a steady  mean-flow  component  and  an  unsteady  fluctuating  component.  Thus, 
for  a general  property  $,  we  have 

4>  - ♦(x.y)  + 4’(x,y,z,t)  , (1) 
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for  a two-dimensional  boundary  layer.  When  expressions  of  this  form  are  sub- 
stituted into  the  equations  of  motion,  the  barred  quantities  lead  to  the 
standard  boundary  lay*r  equations  (after  an  appropriate  ordering  of  terms 
has  been  made),  whereas  the  primed  terms  lead  to  the  stability  equations. 

A complete  derivation  of  these  equations  is  given  in  standard  textbooks  for 
the  incompressible  case;  the  most  complete  drrivation  for  the  compressible 
case  is  given  by  Mack  (1969) . 

The  stability  equations  are  traditionally  solved  by  Laplace  transform 
techniques  in  time,  and  Fourier  transform  techniques  in  the  streamwise 
direction.  For  example,  the  perturbation,  can  be  expressed  as 

♦’(x,y,t)  - $(y)  expQi(ax-wt)  "J  • (2) 

Upon  substitution  of  expressions  of  this  form  into  the  linear  stability 
equations,  we  obtain  a set  of  coupled,  ordinary  differential  equations  in 
the  cross-stream  variable,  y.  The  resulting  set  of  ordinary  differential 
equations  can  then  be  solved  by  standard  numerical  techniques.  Since  these 
ordinary  differential  equations  are  homogeneous,  and  since  the  corresponding 
boundary  conditions  are  likewise  homogeneous  (eee  Mack,  1969,  for  details), 
the  only  solution  which  can  be  obtained  (other  than  a trivial  solution)  is 
sn  eigenvalue  solution.  If  the  frequency,  u,  is  taken  as  purely  real,  while 
the  wave  number,  a,  is  allowed  to  be  complex,  the  spatial  amplification  rate 
of  a wave  of  a particular  frequency  and  Reynolds  number  can  be  obtained  as  the 
eigenvalue  of  the  system.  If  the  amplification  rate  of  a single  frequency 
disturbance  is  computed  at  each  of  a number  of  streamwise  locations  in  the 
boundary  layer  (with  the  boundary  layer  being  assumed  to  be  locally  parallel, 
but  of  a different  thickness  at  each  location),  the  resulting  amplitude  of 
this  wave,  in  terms  of  its  amplitude  at  some  reference  location,  is 

A(o))/Aq(w)  ■ exp  E/a1(x,e,)dx:  , (3) 

where  a dependence  on  the  frequency  is  assumed. 

The  adaptation  of  linear  stability  theory  to  the  problem  of  predicting 
the  significance  of  various  parameters  on  transition  is  achieved  by  tabulating 
the  dependence  of  the  amplification  rate,  -a^,  on  these  parameters  (by  solu- 
tions of  the  linear  stability  equations).  It  should  be  noted  that,  at  this 
Btage,  the  analysis  is  still  completely  analytical;  no  empirical  information 
has  been  Included.  The  response  of  disturbances  of  various  frequencies,  and 
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thelr  net  growth  or  decay  can  be  obtained  directly  from  the  equations  of 
motion.  This  represents  one  of  the  most  attractive  features  of  the  linear 
stability  approach  to  transition;  it  includes  the  physics  of  the  pre- 
transitional  boundary  layer  (so  long  as  disturbances  are  small),  and  it  1b 
simple  enough  to  be  tractable.  Methods  for  applying  these  stability  results 
to  the  actual  prediction  of  transition  are  given  below. 

2.4  Transition  Criteria  for  Use  with  Linear  Stability  Theory 

In  order  to  link  the  linear  stability  theory  with  transition,  some 
empirical  transition  criterion  is  necessary.  If  the  non-linear  effects  were  j.j 

Included  in  the  stability  calculations,  such  an  empirical  criterion  could, 
in  principle,  be  bypassed,  although  it  would  probably  still  be  practical  to 
Include  it.  One  of  the  most  widely  used  techniques  for  predicting  boundary 
layer  transition  ia  the  e9  method  suggested  by  Smith  and  Gamberonl  (1956), 
and  used  in  our  previous  work;  Markle,  Kubota  and  Ko  (1974).  In  the  e9  method, 
the  growth  of  disturbances  is  computed  by  linear  stability  theory  until  an  em- 
pirically defined  leval  of  amplification  has  bean  exceeded,  at  which  time, 
transition  is  said  to  occur.  Although  an  amplification  of  a9  is  reasonable 
in  incompressible  flows,  it  is  considerably  too  large  for  compressible  boundary 
layers,  where  transition  can  frequently  occur  when  the  disturbances,  as  pre- 
dicted by  linear  stability  theory,  have  grown  by  only  a factor  of  e2  (Markle, 

Kubota  and  Ko,  1974),  A potentially  useful  empirical  transition  criterion 
for  compressible  boundary  layers  was  suggested  by  Mack  (1975).  In  this 
approach,  the  initial  disturbance  levels  in  the  boundary  layer  are  related  to 
the  free  stream  disturbances  by  means  of  an  empirical  formula.  Than  the  growth 
of  disturbances  is  computed  by  linear  stability  theory,  with  transition  being 
predicted  when  the  amplitude  of  the  disturbance  reaches  a particular  leval. 

With  this  approach.  Mack  was  able  to  use  a single  transition  criterion  for  all 
supersonic  boundary  layers,  but  he  still  was  forced  to  use  a different  relation 
for  subsonic  boundary  layers.  We  have  not  as  yet  used  this  criterion  in  our 
nosetip  transition  predictions,  but  plan  to  test  it  in  the  future. 


Flow  Research  Report  No.  60 
July,  1975 


3.  PREDICTION  OF  THE  EFFECTS  OF  VARIOUS  PARAMETERS  ON  TRANSITION 

This  section  presents  some  transition  predictions  which  have  been  obtained 
by  applying  the  results  of  linear  stability  theory.  In  order  to  assess  their 
accuracy,  these  predictions  are  compared  with  appropriate  wind-tunnel  measure- 
ments. The  general  objective  of  these  comparisons  is  to  evaluate  the  capability 
of  the  stability  approach  to  predict  the  effects  of  various  parameters  on  the 
location  of  transition. 

3 . 1 The  Effect  of  Reynold;;  Itober-uar-Inch 

Measurements  of  the  location  of  transition  on  a model  in  a supersonic  wind- 
tunnel  have  consistently  shown  that  transition  depends  on  the  Reynolds  number- 
par-inch  of  the  free-stream  flow.  Similarly,  transition  measurements  of  the 
sans  modal  in  different  wind-tunnels  has  shown  substantial  variations.  These 
variations  in  the  transition  location  have  been  frequently  attributed  to  tunnel- 
to-tunnel  differences  in  the  free-stream  disturbance  levels,  or  on  variations 
In  the  disturbance  levels  within  a single  tunnel  with  Reynolds  number-par-inch. 
Xn  fact,  Pate  and  Schueler  (1969)  and  Pate  (1971)  correlated  a wide  range  of 
data  solely  on  the  basis  of  the  acoustic  environment  in  the  wind-tunnel.  Their 
results  suggest  that  the  transition  location  on  a flat  plate  or  a cone  is  inde- 
pendent of  Mach  number,  and  depends  only  on  the  disturbance  environment. 

In  order  to  predict  the  observed  Reynolds  number-per-inch  variations  in 
the  location  of  transition,  the  linear  stability  results  must  be  coupled  with 
a transition  criterion  which  depends  upon  the  initial  disturbance  level  in  the 
boundary  layer. 

The  en  transition  criterion  which  has  been  developed  by  Smith  and  Gam- 
beronl  (1956)  indicates  that  transition  will  occur  when  disturbances  within 
the  boundary  layer  have  grown  by  a specific  (empirically  determined)  factor. 
However,  this  link  between  the  total  amplification  of  disturbances  and  transi- 
tion is  based  upon  an  implicit  assumption  that  the  initial  disturbance  levels 
are  similar  from  case  to  case.  In  situations  in  which  the  initial  disturbance 
level  changes,  the  e11  criterion  (which  is  based  on  totui  amplification)  should 
be  replaced  by  an  amplitude  criterion  which  assumes  that  transition  will  begin 
when  the  disturbance  level  inside  the  boundary  layer  reaches  a particular 
threshold  value.  The  use  of  such  an  amplitude  criterion,  which  explicitly 
accounts  for  the  initial  disturbance  level,  would  lead  to  earlier  transition 
predictions  in  "noisy"  wind  tunnels  than  in  "quiet"  wind  tunnels.  By  inference, 
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the corresponding  amplification  ratios,  A/Aq,  at  the  predicted  transition 
location  will  be  smaller  for  the  noisy  wind-tunnel  case  than  for  the  quiet 
wind-tunnel  case.  We  believe  that  this  interpretation  of  the  en  method  is 
in  the  same  spirit  that  Smith  originally  intended  the  method  to  be  used. 

Before  discussing  the  transition  predictions,  it  is  worthwhile  to  re- 
view the  general  character  of  wind-tunnel  diaturbences. 

It  is  well-known  that  small  disturbances  in  compressible  flows  can  be 
grouped  into  one  of  three  types,  namely,  vorticity  fluctuations,  pressure 
fluctuations,  and  entropy  fluctuations.  Of  these  three  types  of  distur- 
bances, the  vorticity  fluctuations  in  the  wind-tunnel  frea-straam  flow 
dominate  the  affect  of  free-straam  disturbances  on  transition  at  low  Mach 
numbers,  whereas  the  pressure,  or  acoustic,  fluctuations  dominate  the  transi- 
tion beha/ior  at  higher  (above  Mach  2.0)  wind-tunnal  speeds.  Xt  appears  that 
entropy  fluctuations  are  of  minor  importance  except  perhaps  in  the  case  of 
heated-inlet  wind-tunnel  testing. 

Experimental  results  of  the  effect  of  frea-straam  turbulence  on  transi- 
tion Reynolds  numbers  wore  reported  by  Schubauar  and  Klebanoff  (1952)  and  a 
composite  curve  was  given  by  Dryden  (1954).  The  combined  effects  of  turbulence 
and  acoustic  waves  on  the  location  of  transition  has  bean  measured  by  Spangler 
and  Wells  (1969).  The  results  of  these  experiments  are  summarised  in  Fig.  1, 
and  show  that  transition  can  behave  very  differently  under  the  Influence  of 
different  types  of  frea-stream  disturbances. 

Mosetip  transition  tests  are  generally  run  at  high  supersonic  wind-tunnel 
Mach  numbers  where  the  effects  of  acoustic  disturbances  are  generally  more 
significant  than  are  the  effects  of  free-stream  turbulence.  Laufer  (1964) 
has  shown  that  the  acoustic  energy  in  the  free-stream  is  approximately  pro- 
portional to  the  free-stream  dynamic  pressure.  This  acoustic  energy,  which  is 
radiated  by  the  boundary  layers  on  the  wind-tunnel  walls,  also  scales  directly 
with  the  skin  friction.  Thus,  as  the  Reynolds  number-per-inch  of  the  flow  la 
increased,  the  Reynolds  number  of  the  wall  boundary  layer  Increases  and 
(assuming  the  Reynolds  number  is  sufficiently  high  that  the  wall  boundary  layer 
is  turbulent),  the  skin  friction  coefficient  decreases.  As  a result,  the 
acoustical  disturbance  level  in  the  boundary  layer  decreases  as  the  Re/in.  is 
increased.  This  suggests  that  the  transition  Reynolds  number  would  increase 
with  the  tunnel  Re/in.  (Specifically,  if  transition  occurs  when  the  disturbance 
amplitude  roaches  a critical  level,  and  if  the  initial  level  is  lowered,  then 
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more  amplification  la  required  and  transition  moves  to  higher  Reynolds  numbers.) 
The  experimental  results  show  that  thlu  la  the  proper  qualitative  trend. 

An  experimental  measure  of  the  changes  in  the  transition  Reynolds  number 
with  Reynolds  number-par-inch  is  shown  in  Fig.  2.  In  these  experiments,  which 
were  performed  by  Pate  and  Schueler  (1969),  the  free-stream  acoustic  disturbance 
level  was  measured,  and  we  have  used  the  acoustic  disturbance  level  aB  the 
abscissa  in  Fig.  2.  As  can  be  seen,  the  experimental  results  indicate  lower 
transition  Reynolds  numbers  as  the  acoustic  disturbance  increases  In  energy. 

The  corresponding  transition  predictions,  which  have  been  baaed  on  thfc  results 
of  stability  theory,  and  a maximum  amplitude  transition  Criterion,  are  also 
shown  in  the  figure.  (For  these  results,  the  maximum  amplitude  criterion  waa 
selected  to  give  the  proper  transition  Reynolds  number  at  one  point.  Thus,  tha 
level  of  the  predictions  has  been  adjusted  to  fit  the  data ; Only  the  rata  of 
change  of  the  predictions  with  the  acoustic  environment  is  meaningful.)  For 
both  tha  Mm  ■ 3.0  and  Mn  « 5.0  casaa,  tha  predictions  indicate  the  correct 
qualitative  trends,  but  undsrpredict  the  effect.  This  discrepancy  could  be 
attributed  to  a number  of  factors.  For  example,  the  reported  transition 
locations  ere  end-of-transitlon  locations.  Stability  theory  should  not  be 
used  for  and-of-transltlon  predictions;  only  for  beginning-of -transition  pre- 
dictions. Nevertheless,  we  believe  that  the  under prediction  of  the  effect  of 
free-stream  disturbances  is  more  fundamental  than  this.  Stability  theory  in- 
dicates that  transition  is  very  dependent  on  the  frequency  of  the  external 
disturbance  environment  as  well  as  its  level.  An  estimate  of  the  frequency 
spectrum  of  the  disturbances  (which  is  not  recorded)  should  be  included  in  the 
initial  disturbance  levels,  AQ.  Such  an  estimate  could  bo  obtained  from  the 
spectral  measurements  of  Laufer  (1964).  However,  Mack  (1974),  using  other  data 
sourcaa,  has  reported  that  tha  Inclusion  of  the  frequency  dependence  is  still 
not  enough  to  predict  the  observed  variations  in  Re/in.  In  his  latest  work, 
Mack  (1975)  has  shown  that  tha  disturbance  level  and  its  frequency  distribu- 
tion must  be  combined  with  a strongly  non-linear  relationship  between  free- 
stream  disturbances  and  boundary-layer  disturbances  in  order  to  successfully 
predict  the  observed  Ra/in.  effects.  This  relationship  represents  the  recepti- 
vity of  the  boundary  layer  to  disturbances  from  the  free-stream,  and  suggests 
that  they  can  be  damped  or  amplified  as  they  interact  with  the  boundary  layer. 
Similar  conclusions  can  be  reached  by  solving  the  "forced"  problem  of  stability 
theory  (Mack,  1971). 
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3 . 2 Extension  of  Turbulent  Sublayer  Model  to  Compressible  Flow 

In  a previous  report  (Merkle,  Kubota,  and  Ko,  1974),  we  developed  a 
"turbulent  sublayer"  model  to  describe  the  effects  of  distributed  surface 
roughness  on  transition.  This  turbulent  sublayer  model  is  baaed  on  the 
experimental  observations  of  Klebanoff  and  Tidstrom  (1972)  which  indicated 
that  aurface  roughness  affects  transition  in  two  distinct  manners;  (1)  rough- 
ness modifies  the  mean  velocity  profiles,  and  (2)  roughness  increases  the 
initial  disturbance  level  inside  tua  boundary  layer.  The  observed  modifica- 
tions in  the  mean  flow  profiles  decreased  the  critical  Reynolds  number  below 
its  smooth-wall  value,  and  Increased  the  rate  of  amplification  of  disturbances 
in  the  (enlarged)  unstable  region.  When  coupled  with  the  higher  initial  dis- 
turbance levels,  this  increased  amplification  rata  lad  to  substantial  reductions 
In  the  transition  Reynolds  number. 

The  turbulent  sublayer  modal  visualises  the  distortions  In  the  mean  flow 
profiles  as  being  the  result  of  an  anhancod  momentum  transfer  near  the  wall 
which  is  caused  by  an  unsteadiness  in  the  viscous  flow  over  the  individual 
roughness  elements.  The  quantitative  effects  of  this  enhanced  momentum  transfer 
are  Included  in  the  computation  of  the  mean  flow  profiles  by  means  of  an  eddy 
vlucoslty  formulation. 

In  our  previous  report,  this  model  was  used  to  predict  the  location  of 
transition  for  a series  of  incompressible  boundary  layers,  including  both 
favorable  and  unfavorable  pressure  gradient  cases.  Before  being  applied  to 
nosetip  boundary  layers,  the  turbulent  sublayer  model  must  be  extended  to 
compressible  flow.  This  extension  has  been  accomplished  by  introducing  an 
eddy  dlffusivity  to  represent  the  enhanced  heat  transfer  in  the  region  near 
the  rough  surface.  This  eddy  dlffusivity  was  then  defined  In  terms  of  the 
momentum  eddy  viscosity  and  a turbulent  Prandtl  number  (which  was  taken  as 
unity).  An  outline  of  this  development  follows. 
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The  equations  for  a compressible,  laminar  boundary  layer  are  (Hayes  and 
Probstein,  1959), 

& 3 

— pur  + -r—  pvr  ■ 0 (4) 

3x  p o 3y  K o 


3u  . 3u  3p  . 3 / 3u  \ 
pu  3x  + p By  “ " 3x  + 3y  3y/ 


„ . A 3H  ^ kT  3H  . / kT  | 3uz/2 

P 3x  + p 3y  3y  Cp  3y  + y*T  ~ C 3y  ’ 

where  we  have  represented  the  effective  viscosity  and  thermal  conductivity  as 


PT  - v + pe 


k^  ■ k + p^/Pr 


respectively.  Then,  by  defining  an  effective  Prandtl  number, 


WTCp 

PrT  " ~T^  " ?r 


1 + th 

i +^f7v. 


we  can  re-write  the  energy  equation  in  its  familiar  laminar  form; 


3H  , . 

PU  3x  + pV 


3H  . _L  [ PT  3H  n / 1 \ 3u2/2*1 

3y  3y  [PrT  3y  PT  ^ PrT } 3y  J ' 


These  equations  can  bs  reduced  to  similarity  form  by  introducing  the 
Levy-Lees  transformation  (Hayes  and  Probstein,  1959), 


5(x)  “ / p u u r2dx 
j e e e 


u r r 

n(x,y)  - / pdy 
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and  defining  the  new  independent  variables, 

F'  « u/u 


G - H/H  . 
e 


The  resulting  similar  form  of  the  equations  is, 


£cF"J  + FF"  + 0(G-F,5i)  - 0 

fe]  +FG'+7T^[(1-^)  H 


where 


For  our  Incompressible  results,  the  eddy  viscosity  in  the  turbulent  sub- 
layer was  expressed  algebraically  as, 


Vkkjl  - , 


where  K^,  A , and  3j  are  constants.  For  the  compressible  variable  property 
case,  we  retain  this  same  expression,  but  after  non-dimensionalizlng  as  re- 
quired in  £q.  (17),  it  differs  from  the  corresponding  incompressible  formula 
by  a viscosity  ratio.  Thus,  for  compressible  flow,  we  have, 


e vk 

r Kt  v Rfk 


Flow  Research  Report  No.  60 
July,  1975 


-16- 


For  the  compressible  flow  calculations  which  have  been  made  to  date,  the 
constants  in  Eq.  (19)  have  been  kept  at  the  same  values  as  were  used  in  the 
incompressible  solutions.  The  only  additional  constant  which  must  be  evaluated 
for  the  compressible  case  is  the  "turbulent"  Frandtl  number,  PrT,  which  was  set 
equal  to  unity.  To  summarise,  the  constants  for  the  compressible  case  are. 


K - 0.1  9,  “ 1.0 

c * 

A+  - 40  PrT  - 1.0 


(20) 


Finally,  the  Sutherland  viscosity  law  was  used  to  specify  the  viscosity-tempera- 
ture relation. 

The  results  of  some  compressible,  rough-wall  solutions  are  given  in 
Section  4.2. 

3.3  Effects  of  Ablation  on  Transition 

The  effects  of  ablation  on  boundary  layer  transition  have  also  been 
estimated  from  linear  stability  results,  and  ere  Bhown  in  Figs.  3 and  4. 

These  predictions,  which  are  far  an  incompressible  boundary  layer,  are  based 
on  an  approximation  of  the  stability  characteristics  of  tha  blown  boundary 
layer  which  is  analogous  to  the  one  which  was  used  to  predict  the  effects  of 
roughness  on  transition.  The  mean  velocity  profile  for  a typical  blown 
boundary  layer  Is  shown  in  Fig.  3.  This  profile  la  compared  with  two  other 
"equivalent"  profiles.  In  the  first  case,  the  Falkner-Skan  profile,  which 
has  the  same  shape  factor,  H,(«5*/0)  ar.  the  blown  velocity  profile,  is  shown. 

For  the  second  case,  the  Falkner-Skan  profile,  which  has  the  same  shear  stress 
at  the  wall,  tw,  is  shown.  In  order  to  obtain  transition  predictions  without 
the  necessity  of  computing  the  stability  properties  of  the  blown  profiles, 
their  stability  characteristics  were  assumed  to  be  approximated,  in  turn,  by 
each  of  these  two  "equivalent"  profiles.  This  resulted  in  two  approximations 
for  the  transition  Reynolds  number  at  each  blowing  rate.  These  results  are 
shown  in  Fig.  4 in  terms  of  the  predicted  transition  Reynolds  number  as  a 
function  of  the  non-dimensional  blowing  rate,  B’,  for  each  of  three  free-stream 
pressure  gi  dients,  namely,  0-0,  0 - 0.2,  and  0-0.5.  The  0-0  case  cor- 
responds to  a flat  plate,  while  the  0-0.5  case  corresponds  to  the  flow  past 
an  axisymmetric  stagnation  point.  The  solid  lines  represent  the  predicted 
transition  Reynolds  numbers  based  on  the  equivalent  assumption.  As  can  be 
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seen,  both  approximations  give  essentially  the  same  predictions.  These  pre- 
dictions indicate  that  blowing  strongly  decreases  the  transition  Reynolds 
number  for  all  three  pressure  gradient  cases.  It  should  be  noted  that  a 
typical  ablation  rate  for  a re-entry  nosetip  corresponds  to  a blowing  rate 
of  about  B'  *0.5  (or  less).  The  predictions  on  this  curve  extend  to  much 
higher  levels  of  wall  mass-flux  than  are  realistic  in  order  to  indicate  the 
trends,  and  in  order  to  bring  the  transition  Reynolds  number  for  the  highly 
favorable  pressure  gradient  case  down  to  the  level  at  which  stability  results 
were  available.  These  predictions  Indicate  that  the  levels  of  ablation  which 
are  encountered  on  a re-entry  nosetip  are  enough  to  reduce  the  transition 
Reynolds  number  by  about  a factor  of  two  from  the  smooth-wall,  non-blown  case. 
However,  as  noted  below,  blowing  from  a rough  surface  may  give,  qualitatively, 
the  opposite  effect  (at  least  for  small  values  of  blowing).  Finally,  it  should 
be  noted  that  the  flat  plate  and  axlsymmetric  stagnation  point  cases  were  en- 
puted  for  a Prandtl  number,  0.7,  and  for  a viscosity-temperature  relation  of 
the  form,  u*TW,  where  u - 0.7.  The  0-0.2  case  wss  for  a Prandtl  number  of 
unity  in  which  the  viscosity  was  taken  to  be  proportional  to  the  temperature. 

3.4  Comparisons  between  Stability  Theory  Predictions  and  the  Aerotherm 

Correlation 

At  this  point,  is  useful  to  compare  the  general  form  of  the  linear 
stability  predictions  with  the  transition  correlation  which  was  developed 
from  the  PAMT  series  data  by  Anderson  (1973).  His  correlation  related  the 
momentum  thickness  Reynolds  number  at  transition  to  the  surface  roughness, 
and  the  wall-to-free-stream  temperature  ratio  as. 


This  correlation  is  for  wind-tunnel  data  taken  on  non-ablating  models,  but  it 
has  since  been  modified  (Anderson,  1975)  by  means  of  analytical  arguments  so 
that  it  will  include  the  effects  of  ablation.  Based  upon  an  assumption  that 
the  "disturbance  level"  which  exists  inside  the  pre-transltional  boundary 
layer  is  proportional  to  the  kinetic  energy  of  the  flow  at  the  top  of  the 
roughness  element,  Anderson  extended  the  relation  to  the  form 
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(22) 


where  the  temperature  ratio,  has  been  replaced  by  the  density  rqtio, 

p /p  . For  a perfect  gas,  these  are,  of  course,  identical,  but  for  the  non- 
perfect  gee  which  exists  near  the  surface  of  a re-entry  vehicle,  they  can  be 
considerably'  different.  The  argument  that  the  disturbance  level  depends  on 
the  kinetic  energy  ratio  suggests  that  it  is  the  density  ratio,  not  the 
temperature  ratio,  that  is  important  in  the  transition  process. 

In  our  roughness  model,  the  effect  of  the  roughness  enters  in  two  ways: 
through  the  roughness  Reynolds  number,  Re^  - u^k/v,  , and  through  the  length 
ratio,  k/dk,  where  <k  represents  some  (unspecified)  integral  thickness  of  the 
boundary  layer.  That  is  to  aay,  the  roughness  appears  in  two  dimensionless 
forms:  the  characteristic  boundary  layer  thickness,  6^,  and  the  length  scale, 

v^/u^.  The  kinetic  energy  ratio  (which  was  used  by  Anderson)  and  the  rough- 
ness Reynolds  number  are  closely  related.  In  particular,  we  can  express  the 


velocity,  u^,  as 


Uu  " k *7  “ k \ 


(23) 


Then  we  have 


Re,. 


pk"k 


peue 


2 RV  <VV“ 


(24) 


where  the  viscosity  ratio  has  been  replaced  by  a temperature  ratio.  Thus, 
the  two  relations  differ  by  a Reynolds  number  and  a temperature  ratio. 

A qualitative  comparison  of  the  stability  analysis  and  the  transition 
criterion  developed  by  Anderson  is  given  in  Table  I.  The  table  compares  the 
predicted  effects  of  roughness,  blowing,  blowing  in  combination  with  rough- 
ness, and  presBure  gradient.  Both  linear  stability  theory  and  the  Aorotherm 
criterion  predict  that  roughneas  has  a strong  effect  on  the  transition  Reynolds 
number.  The  Aerotherm  criterion  indicates  that  transition  will  never  occur  if 
the  surface  is  perfectly  smooth  (i.e.,  Re^j  + ^aB  fc*0).  Our  stability-based 
roughness  model,  in  turn,  Indicates  that  there  is  a threshold  roughness  height 
below  which  roughneas  has  no  effect  (or  a negligible  effect).  Thus,  in  the 
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limit  as  k+0,  stability  theory  predicts  that  the  transition  Reynolds  number 
will  approach  a (finite)  constant.  Both  of  these  behaviors  suggest  that, 
for  smooth  surfaces,  phenomena  in  addition  to  roughness  control  the  transition 
process.  These  additional  phenomena  have  been  (implicitly)  included  in  the 
stability  approach,  but  have  been  omitted  in  the  Anderson  calculations. 

Because  of  its  limitation  to  finite  values  of  roughness,  the  Aerotherm 
criterion  cannot  be  used  to  predict  the  effect  of  ablation  (blowing)  on  a 
smooth  surface.  The  stability  results,  as  indicated  above,  predict  that 
blowing  has  a strong  de-stabilislng  effect  on  the  boundary  layer;  that  is, 
blowing  moves  transition  rapidly  forward. 

The  combination  of  blowing  in  the  presence  of  roughness  behaveo  quite 
differently  from  blowing  in  conjunction  with  a smooth  surface.  The  Aerotherm 
criterion  indicates  that  blowing  tends  to  diminish  tha  effect  of  the  surface 
roughness  and  causes  transition  to  move  backward  (away  from  the  stagnation 
point).  Definitive  results  from  tha  stability  predictions  have  not  been 
computed  yet,  but  initial  observations,  based  on  qualitative  comparison  of 
the  mean-flow  profiles,  indicate  that  the  stability  results  will  also  predict 
an  initial  increase  in  the  transition  Reynolds  number  for  mild  amounts  of 
blowing.  If  the  blowing  rate  is  increased  to  successively  higher  values,  it 
will  eventually  control  the  transition  location  and  cause  it  to  move  forward 
once  again.  Thus,  the  stability  analysis  indicates  that  there  will  be  a 
roughness-dominated  regime  in  which  increased  ablation  stabilises  the  boundary 
layer  and  causes  the  transition  location  to  move  backward;  there  will  also  be 
a blowing-dominated  regime  in  which  increased  ablation  destabilizes  the 
boundary  layer. 

4.  STABILITY  CALCULATIONS 

Although  we  have  been  able  to  use  existing  stability  calculations  which 
have  been  reported  in  the  literature  to  obtain  some  stability-based  transition 
predictions,  and  to  determine  the  relative  Importance  of  some  parameters  in 
controlling  the  transition  location,  it  is  still  imperative  that  we  generate 
additional  stability  calculations  to  enable  us  to  investigate  the  complete 
spectrum  of  mechanisms  which  enn  affect  the  nosetip  boundary  layer.  In  order 
to  obtain  the  capability  for  computing  these  required  results,  we  have  acquired 
a copy  of  the  linear  stability  code  which  has  been  developed  by  Mack  (1574)  at 
JFL  over  the  past  decade.  Our  version  of  this  program  1b  complete,  and  is  in  a 
"production"  state.  The  program  will  compute  the  stability  properties  of  either 
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two-dimensional  or  three-dimensional  (skewed  waves)  disturbances  in  a boundary 
layer.  It  solves  the  complete,  viscous  stability  equations  for  either  com- 
pressible or  incompressible  flows  and,  in  addition  to  being  capable  of  handling 
two-dimensional  mean  flows,  it  can  also  handle  three-dimensional  boundary  layers. 

Our  version  of  this  computer  program  has  bean  carefully  verified  on  a number- 
to-number  veals  against  appropriate  test  cases  which  were  supplied  by  Mack  to 
ensure  that  all  machlne-to-uachlne  incompatibilities  which  wsre  introduced  by 
the  required  machine  conversion  ware  removed.  Following  this,  the  coda  was 
checked  against  the  incompressible  tabulations  of  Wasson,  Okamura  and  Smith  (1968) 
to  ensure  that  our  version  of  Mack's  code  gave  results  which  were  identical  to 
thosa  which  had  been  obtained  from  other  Independent  codes.  To  ensure  maximum 
verification,  our  "Incompressible"  teat  cases  ware  run  at  a Mach  number  of  0.01 
so  that  the  cods  was  forced  to  use  the  compressible  flow  logic  which  was  to  be 
used  for  the  cases  of  interest  in  the  nosetip  problem.  This  check  also  verified 
that  the  coupling  between  our  mean-flow  boundary  layer  program  (similarity  solu- 
tion) and  the  stability  program  was  proper.  The  results  of  both  xero  and  non- 
Karo  pressure  gradient  calculations  from  our  stability  program  proved  to  be 
identical  to  the  results  which  ware  reported  by  Wssssn  al.  (1968)  (to  the 
accuracy  to  which  we  could  read  thsir  published  curves) . 

Following  the  verification  of  the  code,  some  minor  modifications  were  made 
to  adapt  the  code  to  our  specific  needs.  The  most  significant  modification  was 
to  add  a control  loop  to  the  program  which  would  allow  the  eigenvalues  et  e 
sequence  of  frequencies  on  a constant  Reynolds  number  line  to  be  computed  in 
a single  pass  on  the  machine.  In  this  mode,  the  computer  stores  previous 
eigenvalues  and  uses  them  to  compute  an  accurate  new  guess  for  the  eigenvalue 
at  the  next  higher  (or  lower)  frequency.  Because  of  the  improved  initial 
guesses  which  can  be  obtained,  this  control  loop  not  only  decreases  the  number 
of  person-hours  which  are  required  to  compute  a stability  map,  but  it  also  de- 
creases the  amount  of  computer  time  which  is  required  per  eigenvalue.  A paral- 
lel capability  for  computing  the  eigenvalues  at  a series  of  Reynolds  numbers 
on  a fixed  frequency  line  was  also  Included.  In  addition  to  Incorporating  tills 
new  control  loop,  a few  other  minor  changes  were  made. 

The  stability  computations  which  we  have  made  with  our  code  can  be  divided 
Into  two  specific  categories > a parametric  series  of  calculations  and  some 
actual  nosetip  boundary  layer  calculations.  These  two  types  of  calculations 
are  discussed  separately  below. 
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4.1  Parametric  Stability  Calculations 

The  parametric  stability  calculations  are  being  computed  so  that  a data 
base  of  stability  results  will  be  available  for  use  in  making  rapid,  approxi- 
mate estimates  of  the  effects  of  various  parameters  on  the  transition  location. 
These  parametric  solutions  are  to  be  used  in  a manner  similar  to  that  used  for 
the  Incompressible,  rough-wall  transition  predictions  which  were  reported  pre- 
viously (Merkle,  Kubota,  and  Ko,  1974)  and  for  the  effects  of  ablation  which 
were  described  in  Section  3 of  the  present  Report.  Of  course,  in  the  com- 
pressible boundary  layer,  a two-parameter  family  of  solutions  will  have  to  be 
considered  in  order  to  Include  the  Important  affects  of  hast  transfer.  Although 
sufficient  data  is  still  not  available  to  make  these  approximate  predictions, 
we  anticipate  that  these  two  variables  will  be  the  shape  factor,  H (as  was  used 
for  the  Incompressible  results),  and  the  wall-to-ftaa-atraam  temperature  ratio, 

W 

Since  the  transition  location  generally  occurs  within  the  subsonic  portions 
of  the  nosatip  boundary  layer  (Anderson,  1975),  we  have  restricted  our  para- 
metric stability  tabulations  to  low  Mach  number  regions.  In  particular,  most 
of  the  results  have  been  computed  for  a unlty-Mach  number  boundary  layer.  A 
few  results  have  also  been  tabulated  at  M#  - 0.5,  but  these  results  are  to  be 
used  primarily  to  verify  the  accuracy  of  interpolation  between  the  Mach-one 
results  and  the  incompressible  results  (which  have  bean  reported  by  Was sen  at  at., 
1968).  At  present,  the  adiabatic  wall  case  has  bean  computed  for  the  similarity 
profiles  of  the  Mach-one  boundary  layer  for  a range  of  favorablo  and  non- 
favorable  pressure  gradient  cases.  Tho  results  show  that  even  at  these  low  Mach 
numbers,  substantial  compressibility  effects  appear  in  the  stability  solutions. 
For  example,  the  maximum  amplification  rates  for  a specific  value  of  the  pressure 
gradient  parameter  in  the  Mach-one  boundary  layer  are  about  half  the  correspond- 
ing value  in  the  incompressible  boundary  layer.  Similarly,  the  critical  Reynolde 
numbers  of  the  compressible  boundary  layer  occur  soma  twenty  percent  below  their 
incompressible  counterparts.  A complete  listing  of  the  eigenvalues  which  we  have 
computed  for  the  Mach-one  boundary  layer  are  given  in  Table  II.  Some  summary 
plots  of  the  results  are  shown  in  Figs.  5 through  8. 
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Figure  5 shows  neutral  stability  curves;  that  is,  curves  along  which 
disturbance  amplification  rate  is  equal  to  zero,  for  flows  with  unity  free- 
streaa  Mach  number  and  pressure  gradients  ranging  from  highly  adverse,  P--0.1, 
to  quite  favorable,  6 “ 0.2,  where  6 is  the  usual  Falkner-Skan  presaure  gradient 
parameter.  The  neutral  stability  curve  for  a boundary  layer  near  an  axisyrometric 
stagnation  point  (6  * 0.5)  Is  given  In  Fig.  6b.  As  might  be  expected  in  the 
light  of  incompressible  flow  results,  the  stability  characteristics  are  strongly 
affected  by  pressure  gradient;  the  boundary-layer  becomes  increasingly  stcbie  as 
the  magnitude  of  the  favorable  pressure  gradient  Increases. 

The  effect  of  free-stream  Mach  number  on  boundary  layer  stability  charac- 
teristics is  indicated  in  Figs.  6a  and  6b  which  show  neutral  stability  curves 
for  both  incompressible  and  unity  Mach  number  flows  with  6 - 0.1,0. 2 and  0.5. 

The  critical  Reynolds  numbers  corresponding  to  the  Mach-one  cases  are  signifi- 
cantly smaller  than  those  of  the  corresponding  incompressible  flows. 

Neutral  stability  curves  for  disturbances  in  the  form  of  oblique  waves  skewed 
at  an  angle,  tj>,  from  the  normal  to  the  fres-stream  direction  are  shown  in  Fig.  7. 
For  incompressible  flow,  Squire's  theorem  (see  Rosenhead,  1963)  states  that  the 
neutral  stability  curves  when  plotted  in  terms  of  the  normalized  variables, 

$ - F/cos2^  and  ft  - RcosiJj  where  R - /ke^,  will  coincide  for  all  wave  angles. 

As  a result,  it  follows  that  the  most  unstable  wave,  i.e.,  the  lowest  critical 
Reynolds  number,  is  the  two-dimensional  wave  with  ij)  - 0°.  In  contrast  to  this, 
Fig.  7 shows  that  the  lowest  critical  (non-dlmenslonallzed)  Reynolds  number  in 
the  Mach-one  case  corresponds  to  a 30°  wave;  Squire's  theorem  clearly  does  not 
hold  at  Mach  one. 

Figure  8 shows  the  disturbance  amplification  rate,  a^,  normalized  by  cos 
for  a range  of  ft  at  a fixed  value  of  ft.  Negative  values  of  represent  grow- 
ing and,  hence,  unstable  waves.  As  noted  above,  normalized  curves  for  various 
values  of  $ will  collapse  onto  a single  line  for  an  incompressible  boundary 
layer.  For  the  compressible  case  shown  here,  the  normalized  amplification  rate, 
<*^,  la  greatest  for  a wave  skewed  at  60°  to  the  free-stream  direction,  while 
the  largest  non-normalized  amplification  rate  corresponds  to  the  wave  with 
t|>  ■ 30°.  The  maximum  non-normalized  amplification  rate  for  a wave  at  60°  is 
substantially  smaller  than  the  corresponding  rate  at  either  0°  or  30°.  This 
behavior  is  in  agreement  with  the  results  of  Mack  (1969)  who  showed  that  the 
wave  angle  of  the  most  unstable  disturbance,  in  the  sense  of  maximum  amplifi- 
cation rate  overall  frequencies,  increases  rapidly  with  Mach  number  and  is  in 
the  range  from  55°  to  60°  for  Mach  numbers  greater  than  1.6. 
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1£  these  tabulated  results  are  to  be  used  to  compute  parametric  transition 
predictions,  some  type  of  interpolation  must  be  devised  to  enable  us  to  obtain 
rapid,  accurate  estimates  of  the  eigenvalues  at  intermediate  values  of  the  para- 
meters. If  the  eigenvalue  at  a specific  Reynolds  number  and  frequency  is  de- 
sired for  a Mach-one  boundary  layer,  it  can  be  simply  obtained  by  an  inter- 
polation. However,  for  intermediate  Mach  numbers,  pressure  gradients,  eto.> 
a more  difficult  multiple  interpolation  is  required.  Using  the  Mach  number 
interpolation  as  an  example,  the  interpolation  between  the  Mach-one  results 
and  the  incompressible  results  is  to  be  obtained  from  solutions  which  are 
stored  in  the  non-dimensional  form,  Re/Re__  vs.  where  the  subscript 

CK  CR 

"CR"  refers  to  the  critical  value.  A comparison  of  the  gradient  parameter, 

0 - 0.5,  is  shown  in  Fig.  9.  As  can  be  seen,  these  variablea  appear  to  be 
much  more  desirable  for  use  in  Mach  number  interpolation  than  do  the  original 
(Re,u)  variables  (sag  Fig.  5). 

4.2  Nosatip  Stability  Calculations 

In  addition  to  the  parametric  studies  described  above,  aoma  computations 
of  the  stability  characteristics  of  the  boundary  layer  on  a nosatip  have  also 
been  made.  For  these  initial  computations,  a local  similarity  approximation 
was  used  to  compute  the  mean  flow  profiles  and  to  include  the  surface  rough- 
ness effects.  Additional  computations,  which  will  be  based  on  a finite-dif- 
ference solution  of  the  boundary  layer  aquations,  are  planned.  These  nosatip 
calculations  combine  the  affects  of  all  the  individual  parameters  into  a 
single  computation.  In  particular,  realistic  values  of  roughness,  wall- 
temperature,  pressure  gradient  and  local  Reynolds  number  have  been  used.  In 
order  to  determine  the  relative  effects  of  the  various  parameters,  some  smooth- 
wall  and  some  adiabatic-wall  calculations  have  also  been  completed. 

The  values  of  the  nosatip  parameters  which  were  used  in  the  calculations 
have  bean  taken  from  a SNAP  computer  output  which  was  supplied  to  us  by  the  1 

Aero therm  Corporation.  The  computed  properties  are  for  a nosatip  at  an  alti-  j 

tude  of  59,000  feet,  a free-stream  Mach  number,  Ma  ■ 22.6,  and  a free-stream 
total  pressure  of  50.7  atmospheres.  The  Bblatlvo  nosetip  radius  was  initially  | 

1.5  inches,  and  for  the  conditions  just  described,  it  remains  nearly  unchanged  j 

in  sise  an'*,  shape  (hemispherical).  Values  of  wall  temperature,  local  total 
pressure,  local  pressure  gradient  and  Mach  number  were  obtained  from  the  computer 
print-out.  Corresponding  Integral  boundary  layer  properties,  including  the 
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momentum  thickness,  6,  and  the  displacement  thickness,  6*,  were  also  supplied 
to  us  by  Aero therm  from  their  BLIMP  boundary  layer  code.  The  roughness  of  the 
nosatlp  surface  was  estimated  to  be  about  0.6  mils. 

A plot  which  shows  the  variation  In  the  Mach  number  at  the  boundary  layer 
edge,  the  wall-to-free-stream  temperature  ratio,  the  static  pressure  distribu- 
tion and  the  Reynolds  number  based  on  momentum  thickness  as  a function  of  dis- 
tance from  the  stagnation  point  is  shown  in  Fig.  10.  The  static  pressure 
variation  in  Fig.  10  has  also  besn  re-plotted  in  Fig.  11  in  terms  of  the  similar 
pressure  gradient  parameter,  0,  as  defined  in  the  Levy-Lees  transformation. , As 
can  be  seen  in  Fig.  11,  the  pressure  gradient  parameter,  0,  is  nearly  constant 

at  the  value  of  0 - 0.5  for  some  distance,  and  then  begins  to  increase  somewhat. 

• • , '1 
For  reference,  the  pressure  distribution  as  computed  by  modified  Newtonian 

theory  is  also  indudad  in  Fig.  11.  Finally,  it  should  be  noted  that  the  Aero- 

therm  transition  criterion  (Anderson,  1975)  indicated  that  transition  would 

occur  at  about  0.7  inches  downstream  of  the  stagnation  point. 

The  mean  flow  profiles  ware  computed  from  the  Levy-Lees  transformed  equa- 
tions of  motions  (Hayes  and  Frobstein,  1959).  Effects  of  roughness  ware  Included 
by  means  of  the  compressible  version  of  the  turbulent  sublayer  model  which  was 
described  in  Section  3.2.  Then,  the  similar  profiles  were  computed  end  scaled 
to  give  the  proper  momentum  thickness  Reynolds  number  (as  determined  from  the 
BLIMP  calculation). 

At  the  present  time,  stability  characteristics  have  been  computed  at  two 
points  on  the  boundary.  One  point  is  at  an  arc-length  distance  of  S ■ 0.1  in., 
which  is  near  the  stagnation  point.  The  other  point  is  at  a location  of  S ■ 0.7  in 
which  is  near  the  indicated  transition  location.  At  the  first  point,  S - 0.1  in., 
the  important  parameters  of  interest  have  the  following  values t edge  Mach  number, 

Ma  * 0.09;  pressure  gradient,  0 “ 0.50$  and  wall  temperature  ratio,  Tw/Tg  - 0.52. 

A plot  of  the  mean  velocity  distribution  for  the  smooth  wall  case  and  a 0.3  mil 
roughness  case  is  shown  in  Fig.  12.  The  computed  Reynolds  number  based  on  the 

i 

boundary  layer  thickness,  6 (defined  bb  the  location  where  u/u  « 0.999),  is 
R«g  ■ 155.  The  roughness  Reynolds  number,  Re^,  for  this  case  is  Re^  ■ 20.  Be- 
cause of  this  small  value,  the  effect  of  the  roughness  on  the  stability 
properties  at  this  location  wsb  quite  small.  As  a result,  our  analysis  showed 
that  all  disturbances  which  were  introduced  into  the  boundary  layer  were  damped, 
even  in  the  presence  of  roughness.  As  an  indication  of  the  margin  of  stability, 
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an  approximate  value  of  the  critical  Reynolds  number  was  determined.  For  the 
u present  edge  conditions,  but  with  an  adiabatic  smooth  wall,  the  critical  Rey- 

nolds number  is  about  Re^  » 37,000.  Computations  for  a smooth,  cold  wall  were 
not  made,  but  it  Is  known  that  the  smooth-wall  boundary  layer  becomes  more 
stable  with  cooling.  Thus,  it  Is  expected  that  the  smooth  cold-wall  critical 
Reynolds  number  will  be  above  37,000.  Similar  stability  calculations  for  a 
rough,  adiabatic  wall  case  Indicated  a critical  Reynolds  number  of  about  13,000 
based  on  boundary  layer  thickness.  It  is  expected  that  the  roughness  is  suffi- 
ciently large  to  lower  this  value  somewhat  as  the  wall  is  cooled  to  the  indi- 

i > cated  T /T  “ 0.52  conditions,  but,  again,  a verification  of  the  temperature 

w e 

effect  was  not  made.  It  should  ba  noted  that  these  estimates  of  the  critical 
Reynolds  number  are  only  approximate  and  are  intended  to  demonstrate  the  ex- 
cessive margin  of  stability  which  exists  at  thiB  near-stagnation  point  location, 
i Besides  the  stability  of  the  boundary  layer,  the  other  important  point  to  note 

Is  that  the  roughness,  even  at  these  low  Reynolds  numbers,  decreases  the  critical 
Reynolds  number  by  about  a factor  of  three. 

Additional  computations  of  the  stability  characteristics  of  the  boundary 
layer  at  a region  very  near  the  transition  location  (S  » 0.7  in.)  are  currently 
in  progress,  but  have  not  been  completed  as  yet.  The  boundary  layer  edge  con- 
dition at  this  location  ares  M « 0.66,  3 **  0.63,  and  Re.  » 595.  Because 

• 0 

of  the  higher  local  Reynolds  number  at  this  point  in  the  boundary  layer,  the 
G roughness  Reynolds  number,  Re^»  is  also  larger  and  is  expected  to  have  a more 

dominant  affect.  The  computed  value  of  Re^  is  about  150.  For  the  adiabatic 
smooth-wall  case,  the  critical  Reynolds  number  is  around  40,000,  which  is 
somewhat  different  than  the  S - 0.1  case  described  above  because  of  tha  differs 
G ence  in  the  Mach  number  and  pressure  gradient.  Again,  cold-wall  values  are 

expected  to  ba  somewhat  larger.  The  corresponding  stability  maps  for  the  rough- 
wall  cases  have  not  bean  completed  yet.  Thus,  at  the  present  time,  it  is  un- 
certain as  to  whether  or  not  the  combined  effects  of  roughness  and  cooling  will 
decrease  the  critical  Reynolds  number  to  a value  below  the  expected  transition 
Reynolds  number.  The  above  figures  indicate  that  roughness  has  to  decrease  the 
critical  Reynolds  number  by  a factor  of  100  in  order  for  stability  theory  to  be 
able  to  predict  transition  at  the  correct  position.  Comparisons  between  rough 
and  smooth  surfaced  nosetips  suggest  that  roughness  can  easily  have  an  effect 
which  is  this  Btrong  (see  discussion  in  the  next  Section). 


Flow  Research  Report  No.  60 
July,  1975 

-26- 

Without  the  aid  of  the  final  rough-wall  nosetlp  calculations,  it  is  dif- 
ficult to  speculate  on  the  outcome  of  the  predictions;  however,  a few  philo- 
sophical points  are  in  order.  If  the  current  roughness  model  does  not  lower 
the  critical  Reynolds  number  sufficiently  to  enable  transition  to  be  predicted 
at  the  Reynolds  numbers  which  have  been  observed  experimentally,  this  can  be 
easily  rectified  by  an  appropriate  change  in  the  values  of  the  constants  in 
the  turbulent  sublayer  model.  That  is  to  say,  we  could  choose  to  evaluate 
the  constants  in  the  roughness  model  on  the  basis  of  nosetlp  experiments,  not 
flat  plate  data  as  has  been  originally  used.  We  prefer,  however,  to  view  the 
roughness  constants  as  being  fixed.  Consequently,  if  this  conjectured  inabil- 
ity to  predict  nosetlp  transition  does  occur,  we  prefer  to  return  to  the 
equations  of  motion  end  carefully  re-evaluate  the  terms  which  have  been  neg- 
lected in  an  attempt  to  find  specific  effects  which  can  substantially  effect 
the  stability  characteristics  on  the  nosetlp. 

5.  ADDITIONAL  TRANSITION  MECHANISMS 

The  success  (or  failure)  of  any  theoretical  procedure  for  the  prediction 
of  transition  on  a rough  nosetlp  is  tempered  by  the  fact  that  in  including  the 
effects  of  surface  roughness,  some  empirically  determined  constants  must  be 
introduced  in  addition  to  those  which  are  required  for  the  prediction  of  transi- 
tion on  smooth  wells.  Thus,  a "proper"  choice  of  the  constants  will  always 
allow  the  correct  transition  location  to  be  predicted,  or  at  least,  the  theory 
can  be  "tuned"  by  evaluating  the  constants  from  one  series  of  experimental  date 
end  then  using  it  to  predict  other  similar  series.  Although  such  a procedure 
can  and  must  be  used  to  some  extent  to  describe  a process  as  complex  as  the 

nosetlp  transition  phenomena,  it  tends  to  obscure  the  question  of  whether  or 

* 

not  the  proper  mechanisms  are  being  simulated.  One  way  to  verify  more  fully 
the  mechanisms  which  lead  to  transition  is  to  use  the  theory  to  predict  nosetlp 
transition  in  situations  in  which  surface  roughness  is  not  a factor,  i.o.>  to 
consider  smooth-wall  noaafcip  transition.  Then,  if  proper  transition  mechanisms 
can  be  predicted  for  the  smooth-walled  case  (for  which  a theoretical  model 
would  require  a minimal  amount  of  empiricism),  the  rough-wall  predictions  can 
be  accepted  with  more  confidence.  Unfortunately,  this  strategy  loads  to  other 
complications  since  new  mechanisms  take  over  when  roughness  is  eliminated. 
Nevertheless,  it  is  useful  to  consider  the  transition  locations  which  have 
been  observed  on  smooth-walled  nose  tips. 
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5 . 1 Experimental  ''Smooth-Wall'1  Nose  tip  Transition 

Since  it  is  impossible  to  obtain  a surface  which  is  "smooth"  in  the 
mathematical  sense,  we  will  use  the  term  "smooth"  to  refer  to  a surface 
whose  roughness  has  a negligible  effect  on  the  location  of  transition. 
Unfortunately,  the  experimental  data  do  not  show  any  clear-cut  level  below 
which  the  wall-roughness  ceases  to  affect  transition.  However,  a number 
of  transition  tests  on  highly  polished  surfaces  have  been  reported.  Garland 
and  Chauvin  (1957)  reported  transition  Reynolds  numbers  (based  on  momentum 
thickness)  of  between  160  to  600  for  an  8-inch  diameter  sphere-cylinder  whose 
surface  was  polished  to  a finish  of  less  than  25  microinches  of  roughness. 

These  transition  measurements  were  obtained  in  free-flight  experiments  at 
Mach  numbers  between  2.0  and  4.0,  and  at  wall-to-free-atream  temperature 
ratios  of  about  T/T  * 1.0-1. 2.  In  comparison  with  the  25  microinch  sur- 

W © 

face  roughness,  the  laminar  boundary  momentum  thickness  on  the  nosetip  was  on 

the  order  of  one  mil.  Tests  of  a hemisphere-cone  at  similar  flight  conditions 

(Chauvin  and  Spaagle,  1957)  and  with  a similar  25  microinch  roughness  on  the 

nosetip  showed  similar  nosetip  transition  Reynolds  numbers  (Re  * 300  - 400). 

In  view  of  the  small  value  of  the  ratio  of  surface  roughneaB  to  boundary  layer 

thickness  on  these  models,  the  surface  roughness  would  not  be  expected  to  have 

an  important  effect  on  the  transition  location;  however,  the  results  of  a 

second  series  of  teats  with  even  smoother  nosetips  indicates  that  this  is  not 

true.  Buglla  (1957)  and  Hall,  Speegle,  and  Piland  (1957)  tested  nosetips 

whose  surface  roughness  was  about  5 microinches.  Their  results  Indicated  that 

transition  on  the  5 microinch  nosetip  occurred  at  Reynolds  numbers  which  were 

about  a factor  of  four  higher  than  the  corresponding  25  microinch  results.  For 

the  5 microinch  nosetip,  transition  Reynolds  numbers  of  from  800  to  2,200  on  a 

sphere-cylinder  and  from  900  to  1,200  for  the  sphere-cone,  respectively,  were 

observed.  Thus,  it  appears  that  even  for  these  highly  polished  surfaces,  the 

roughness  exerts  a significant  control  over  the  transition  location. 

Additional  evidence  of  the  sensitivity  of  the  location  of  transition  to 

the  surface  finish  on  the  nosetip  was  reported  by  Dunlap  and  Kuethe  (1962). 

They  observed  nosetip  transition  Reynolds  numbers  which  were  around  Re^  - 600 

for  a ratio  of  surface  to  free-stream  temperatures  ranging  from  Tw/Te  * 0.5 

to  T /T  - 1.2,  in  a simulated  hypersonic  environment.  Their  results  were  for 
we 

a "very  highly  polished"  surface  ("probably  around  a few  microinches  of  rough- 
ness"), but  they  noted  that  dust  particles  tended  to  collect  near  the  stagnation 
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point  (within  10  degrees)  and  cause  premature  transition.  Meaningful  data 
could  only  be  obtained  by  frequent  polishing  and  inspection  of  the  nosetip 
for  dust  accumulation,  pitting,  and  frost  due  to  freezing  of  COg.  Thus,  the 
experimental  determination  of  the  transition  location  on  a smooth  nosetip 
may  not  only  be  sensitive  to  how  smooth  the  nosetip  is,  but  also  to  the  care 
with  which  the  experiment  la  conducted.  In  tests  In  which  very  highly  polished 
surfaces  are  used,  the  experimental  conditions  must  be  extremely  closely  moni- 
tored to  ensure  unambiguous  results. 

5 . 2 Additional  Transition  Mechanisms  and  Possible  Extension  of  Stability  Theory 
Stability  theory  indicates  that  the  critical  Reynolds  number  in  axisymmetric 
boundary  layer  near  a stagnation  point  occurs  at  a value  of  Re  ■ 3,300,  which  is 
considerably  higher  than  the  transition  Reynolds  numbers  which  are  quoted  above 
for  highly  polished  surfaces.  Whether  or  not  transition  on  a completely  smooth 
surface  would  occur  at  Reynolds  numbers  which  are  above  this  value  remains  an 
open  question.  Consequently,  it  is  imperative  to  review  the  known  theoretical 
arguments  for  other  potential  transition  mechanisms  in  an  attempt  to  ascertain 
their  effect  on  nosetip  transition.  Thus,  we  approach  the  problem  from  the 
following  viewpoint.  Transition  i3  experimentally  observed  to  occur  on  the 
nosetip  in  regions  where  classical  parallel  flow  stability  theory  (without  the 
aid  of  a roughness  model)  indicates  thst  all  disturbances  are  damped.  If  tran- 
sition, in  fact,  does  occur  as  a result  of  the  growth  of  small  disturbances, 
some  other  physical  phenomena  which  have  been  hitherto  ignored  must  be  included 
in  the  problem  formulation.  Specifically,  some  terms  in  the  stability  equations, 
which  have  been  assumed  to  be  of  higher  order  and  ignored,  may  in  reality  have 
a significant  effect.  One  path  around  this  dilemma  is  through  the  use  of  a 
roughness  model  which  affects  the  mean  velocity  profiles.  However,  this  rough- 
ness analysis  can,  at  most,  be  included  in  an  approximate  way,  and  it  must  in- 
clude empirical  "constants".  To  be  sure,  these  constants  can  easily  be  adjusted 
so  that  instability  and  transition  are  predicted  to  occur  on  the  nosetip,  but 
such  an  arbitrary  adjustment  of  constants  in  order  to  match  the  experimental 
data  has  little  meaning,  especially  if  we  are  interested  in  elucidating  the 
mechanisms  which  lead  to  nosetip  transition.  An  alternative  approach  is  to 
review  carefully  the  ordering  of  terms  in  the  linearization  and  simplification 
of  the  stability  equation  to  determine  which  (If  any)  effects  have  been  inad- 
vertently relegated  to  higher  order  and,  hence,  dropped  from  the  analysis, 
despite  their  first-order  effect  in  the  physical  situation.  Some  of  the  effects 
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which  have  been  neglected  In  our  current  analysis  are: 

1.  Effects  of  streamline  curvature  In  the  approach  to  the  stagnation 
point  which  could  lead  to  Taylor-GBrtler  (longitudinal)  instability; 

2.  Additional  axisymmetric  terms  in  the  equations  of  motion; 

3.  Coupling  between  disturbances  and  the  mean  flow  profiles  (even  in 
regions  where  disturbances  are  damped)  which  act  in  a cumulative  sense 
alter  the  mean-flow  profiles  and,  hence,  alter  the  stability 
characteristics ; 

4.  Unsteadiness  in  the  mean  flow  caused  by  "wandering"  of  the 
stagnation  point; 

5.  Effects  of  non-zero  velocity  gradients  (rotational  flow)  outside  the 
boundary  layer  caused  by  entropy  swallowing  effects. 


A brief  review  of  these  supposedly  "higher  order"  influences  is  given  below. 

Some  evidence  of  the  importance  of  Gttrtler  instabilities,  or  axlsymmetric 
vortex  stretching  effects  was  reported  by  Kuethe  (1957),  who  observed  the  growth  of 


wave-like  disturbances  in  the  boundary  layer  over  a hemisphere  at  local  Reynolds 
numbers  which  were  substantially  below  the  critical  value.  Motivated  by  these 
observations,  he  performed  a simple  analysis  which  Indicated  that  the  vortex 
stretching  effects,  which  occur  as  the  fluid  in  the  boundary  layer  is  swept 
over  the  sphere,  are  responsible  for  the  observed  amplification.  His  model  is, 
however,  essentially  a dimensional  analysis,  and  cannot  by  itself  yield  quanti- 
tative results.  Kuethe'a  analysis  suggests  that  the  stretching  of  a longi- 
tudinal vortex  In  a region  in  which  the  velocity  gradient,  Su/8x,  is  constant 
is  given  by 


v;/v; 


/T77- 


where  represents  the  position  of  a particle  in  the  boundary  layer  at  some 

initial  time,  and  x2  represents  its  position  an  incremental  time  later.  v2 

and  Vj  are  the  corresponding  disturbance  velocities.  A similar  analyBio  for 

the  stretching  of  nu  axlmuthnl  vortex  filament,  eorresnond Lng  to  an  axisym- 

metric  Tollmiun-Schlichting  wave  (which  was  suggested,  but  not  carried  out 

by  Kuethe)  gives,  

vj/v|  - /R?/Rj  , 
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where  R represents  the  radial  coordinate  of  the  sphere.  A quantitative  evalua- 
tion of  these  vortex  stretching  mechanisms  would  require  either  an  analysis  of 
the  complete  Gortler  instability  equations  (for  the  former  case)  or  of  the  axi- 
symmetric  form  of  the  Tollmien-Schlichtlng  equations  (for  the  latter  case).  An 
assessment  of  th*.  feasibility  of  such  an  analysis  is  currently  underway. 

The  inclusion  of  the  coupling  and  the  cumulative  Interaction  between  the 
disturbances  and  the  mean  flow  could  be  easily  incorporated  into  a stability 
analysis  if  only  single  frequency  disturbances  were  considered.  Thu;:.,  if  a 
disturbance/mean  flow  energy  transfer  were  allowed,  the  disturbances  would  feed 
energy  into  the  mean  flow  as  they  decayed  (in  the  stable  region),  and,  as  a 
consequence,  the  mean  flow  profiles  would  become  distorted.  These  distorted 
profiles  could,  in  turn,  generate  an  instability;  however,  the  likelihood  of 
such  a process  being  significant  seems  remote. 

Unsteadiuess  in  the  mean  flow  profiles  has  been  reported  by  Kuethe, 
Willmarth,  and  Crocker  (1960),  who  described  experimental  evidence  of  random 
fluctuations  or  "wandering"  of  the  stagnation  point  about  its  equilibrium 
position  on  hemispherical  wind-tunnel  modelB  in  both  subsonic  and  supersonic 
free-stream  flows.  Similar  observations  have  also  been  reported  by  Peterson 
and  Horton  (1959).  A Gdrtler-like  stagnation  point  instability  (see  above ) 
and  stretching  of  vortex  filament  in  the  diverging  stagnation  point  flow  were 
suggested  as  possible  generation  mechanisms  for  these  mean  flow  fluctuations. 
Such  stagnation  point  "wandering"  would  clearly  lead  to  mean  flow  fluctuations 
which  have  not  been  considered  in  the  present  study. 

The  strongly  curved,  hypersonic  bow  shock  causes  an  inviscid  shear  flow 
layer  to  be  generated  in  the  region  outside  the  boundary  .layer.  In  all  of 
our  calculations  to  date,  this  shear  layer  has  been  Ignored,  and  the  flow  out- 
side the  boundary  layer  has  been  taken  as  being  uniform  and  parallel.  The 
inclusion  of  the  characteristics  of  this  shear  layer  in  the  stability  solutions 
could  be  accomplished  by  simply  extending  the  mean  flow  profiles  which  are  used 
in  the  stability  calculation  beyond  the  "edge"  of  the  boundary  layer.  Again, 
order-of-magnitude  estimates  of  the  effect  of  this  shear  layer  suggest  that  it 
is  not  sufficiently  strong  to  materially  affect  the  stability  characteristics; 
however,  a more  rigorous  analysis  may  be  justifiable. 
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5 . 3 Comparison  with  Turbulent-Modeling  Approaches 

In  conjunction  with  the  abovu  discussion,  it  should  be  pointed  out  that 
several  of  the  items  in  the  list  of  effects  which  are  omitted  in  the  stability 
approach  are  included  in  the  turbulence-modeling  approaches  which  are  being 
used  to  predict  transition  (ees  Section  2.2).'.  Fop  example,  the  turbulence- 
modeling  approaches  do  include  the  energy  exchange  between  the  disturbances 
and  the  mean  flow.  Nevertheless,  it i appears  that  the  energy  loss  from  the 
mean  flow  to  the  disturbances  is  not  a significant  effect  until  after  the 
disturbances  have  grown  quite  large.  That  is  to  aay,  it  appears  that  the 
initial  growth  of  "turbulence",  as  computed  by  these  models,  would  not  be 
affected  If  the  energy  exchange  between  the  mean  flow  and  the  turbulence 
were  neglected.  Thus,  we  believe  that  title  effect  which  1b  included  in  the 
turbulence  models  and  is  ignored  in  the  stability  approach  is  not  Important. 
This  speculation  could  be  checked  for  the  turbulence  models  by  means  of  two 
back-to-back  computations,  the  first  of  which  would  monitor  the  initial  stages 
of  the  growth  of  the  turbulence  in  a standard  computation,  while  the  second 
would  monitor  the  turbulence  growth  in  a computation  with  the  coupling  term 
"unhooked"  from  the  mean  flow,  to  determine  the  downstream  location  at  which 
the  two  turbulence  growth  rates  began  to  deviate  from  each  other. 

A second  difference  between  the  stability  technique  and  the  turbulence 
model  approaches  is  that  the  latter  Include  the  full  axiBymmetric  effects. 

This  could  be  significant,  especially  very  near  the  stagnation  point  where 
the  axlsymmetrlc  terms  dominate.  As  Indicated  above,  an  assessment  of  the 
importance  of  theBe  effects  and  means  of  including  them  in  the  stability 
equations  is  currently  being  undertaken. 

The  longitudinal  (Cdrtler)  vortex  effect  represents  a three-dimensional 
( i.e.y  a variation  with  azimuthal  angle,  <j>)  effect  which  is  ignored  in  both 
approaches.  One  deduction  Is  that  if  it  is  not  Included  in  the  turbulence 
model  approaches,  and  if  they  can  accurately  predict  nosetlp  transition,  then, 
by  implication,  it  is  not  important  in  the  stability  analysis  either.  This 
supposition,  however,  presupposes  that  the  turbulence  models  predict  transition 
at  the  right  position  for  the  right  reason;  in  view  of  the  large  number  of 
empirical  constants  involved,  such  is  probably  not  a justifiable  conclusion. 

We  believe  these  3-1)  effects  could  still  be  important,  perhaps  in  a manner 
which  is  analogous  to  the  observed  spanwise  variations  which  were  observed  in 
planar  boundary  layers,  as  discussed  in  Section  2. 
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Finally,  it  should  be  noted  that  stability  theory  Includes  some  effects 
which  the  turbulence  models  ignore.  One  primary  difference  is  that  stability 
theory  solves  the  complete  time-dependent  equations  for  the  disturbances, 
whereas  the  turbulence  model  approaches  include  only  the  time-averaged  equations. 
Of  course,  the  dominant  difference  is  that  stability  theory  solves  the  actual 
equations  of  motion  while  the  turbulence  models  solve  a set  of  equations  which 
"model"  the  equations  of  motion.  Xn  addition,  these  "models"  have  been  developed 
for  turbulence,  not  for  the  very  different  type  of  unsteadiness  which  character- 
ises the  transitional  boundary  layer. 

Another  item  of  interest  in  the  comparison  between  the  two  techniques  is 
their  treatment  of  the  "production"  and  "dissipation"  of  disturbances.  Xn 
stability  theory,  the  growth  (or  decay)  rate  of  disturbances  is  very  much  de- 
pendent upon  their  frequency  content,  and  the  amplification  rate  for  the  unstable 
frequencies  is  strongly  affected  by  the  existence  of  a point  of  inflection  in  the 
boundary  layer  profile.  Thus,  for  example,  if  we  consider  the  family  of  Falkner- 
Skan  profiles,  the  favorable  pressure  gradient  cases,  which  have  no  point  of 
inflection,  become  increasingly  stable  to  small  disturbances  as  the  pressure 
gradient  parameter  is  Increased.  By  way  of  contrast,  the  unfavorable  pressure 
gradient  cases  have  a point  of  inflection  which  moves  steadily  away  from  the 
wall  as  the  value  of  3 get  more  negative,  and,  correspondingly,  the  amplifica- 
tion gets  increasingly  more  rapid.  Turbulence  models,  however,  generally  have 
a production  term  which  depends  on  the  magnitude  of  the  velocity  gradient. 

This  suggests  that  steep  velocity  gradients  (such  as  occur  in  fully  turbulent 
boundary  layers  or  favorable  pressure  gradient,  laminar  boundary  layers)  in- 
crease the  production  of  turbulence.  Thus,  it  appears  that  a favorable  pressure 
gradient  (positive  value  of  3)  would  lead  to  Increased  production,  whereas  an 
unfavorable  pressure  gradient  would  decrease  the  production  of  turbulence.  Of 
course,  the  net  change  of  turbulence  must  Include  the  dissipative  effects  also, 
and  it  is  the  balance  between  these  two  terms  which  yields  the  net  growth  rate 
of  the  disturbances  in  the  laminar  boundary  layer. 
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6.  SUMMARY  AND  CONCLUSIONS 

i j Linear  stability  theory  has  been  used  as  the  basic  theoretical  foundation 

for  an  analysis  of  the  mechanisms  which  lead  to  boundary  layer  transition  on 
the  noaetip  of  a re-entry  vehicle.  In  this  analysis,  the  effects  of  surface 
roughness,  wall  cooling,  ablation,  pressure  gradient,  and  axisymmetric  geo- 

, metry  have  been  considered  both  individually  and  in  combination  with  each  other. 

The  results  suggest  that  surface  roughness  is  an  important  controlling  parameter 
in  determining  the  location  of  transition,  but  that  the  other  parameters  can 
substantially  alter  the  "effective"  roughness  height. 

I j Some  specific  observations  and  conclusions  which  can  be  made  regarding 

nosetip  transition  includet 

1.  A number  of  physical  processes,  which  have  a profound  Influence  on 
the  location  of  transition,  occur  prior  to  the  "beginning"  of  the  transition 

' .i  region  as  it  is  microscopically  observed.  The  use  of  linear  stability  theory 

in  a technique  for  predicting  transition  is  attractive  because  it  includes 
the  cspsbillty  for  fslthfully  predicting  these  pre-transitlonal  phenomena; 

2.  Noaetip  transition  is  experimentally  observed  to  occur  at  local 
Reynolds  numbers  which  are  substantially  smaller  than  the  critical  value  which 
is  predicted  by  stability  theory.  Consequently,  classical,  parallel-flow 
stability  results,  without  modification,  are  unable  to  explain  these  experi- 
mental transition  locations  (however,  000  Item  3 following); 

3.  The  effects  of  surface  roughness  have  been  Included  in  the  stability 
analyses  of  noaetip  boundary  layers  by  means  of  a "turbulent  sublayer"  model 
which  altera  the  mean  flow  profiles  and  reduces  the  critical  Reynolds  number 
sufficiently  to  allow  the  application  of  stability  theory  in  nosetip-like 
environments.  Some  specific  predictions  of  the  roughness  model  are: 

a.  Experimentally  observed  effects  of  roughness  on  transition 
can  be  predicted. 

b.  The  effect  of  cooling  a rough  surface  tends  to  destabilise 

the  boundary  layer  and  cause  earlier  transition  (because  the  boundary 
layer  thinning  effect  which  is  caused  by  cooling  makes  the  roughness 
more  effective).  By  contrast,  cooling  a smooth  wall  makes  the  bounday 
' ->  > more  stable. 
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c.  The  effects  of  ablation  on  a smooth-wall  boundary  layer  are  to 

^ decrease  the  transition  Reynolds  number.  Preliminary  results  indicate 

that  ablation,  in  combination  with  rough  surfaces,  can,  again,  have  the 

opposite  (i.e. , stabilizing)  effects 

d.  Sufficiently  large  roughnesses  can  essentially  negate  the 

v , strongly  stabilizing  effect  of  highly  favorable  pressure  gradients  such 

as  the  ones  which  occur  in  the  nosetip  region; 

4.  The  effects  of  wind-tunnel  Re/in.  on  the  transition  Reynolds  number 
have  been  investigated,  but  it  appears  that  the  frequency  content  of  the 

, , free-stream  disturbances  must  be  included  before  this  effect  can  be  properly 

predicted; 

5.  An  extensive  tabulation  of  the  stability  characteristics  of  a Mach- 
one  boundary  layer  (for  both  favorable  and  unfavorable  pressure  gradients)  has 
been  compiled  and  is  presented  herein.  These  resultB  show  that,  even  at  unity  Mach 
number,  compressibility  has  a significant  effect  on  the  stability  characteristics; 

6.  Computations  of  the  stability  properties  of  an  actual  nosetip  boundary 
layer  have  indicated  that  regions  near  the  stagnation  point  are  stable  to  all 
disturbances,  even  in  the  presence  of  roughness.  Analyses  further  away  from 
the  nosetip  are  currently  in  progress; 

7.  A number  of  additional  factors  which  could  Influence  the  growth  of 
disturbances  in  the  nosetip  boundary  layer,  but  which  have  until  now  been 

U ignored  in  the  stability  analysis,  have  been  identified  and  are  being  more 

carefully  evaluated. 
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Predicted  effect  of  Reynolds  Haber— per— Inch  on  Transition  In  Hind 
Tunnel  (end  of  transition). 


Flow  Research  Report  No.  60 
July,  1975 


Fig.  3 Comparison  of  Blown  Boundary-Layer  Velocity 
Profile  with  Equivalent  Non-Blown  Profiles. 


Non-Dimensional  Frequency 
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Fig.  6.b  Neutral  Stability  Curves  at  M » 0 and  M 
Favorable  Pressure  Gradient,  6 * 0.5 


NORMALIZED  FREQUENCY*  F/F 


1.2 


iriations  of  Boundary-Layer  Edge  Conditions  around 
a Nosetip. 


Fig.  11  Variation  of  Pressure  Gradient  Parameter  around  a Nosetip 


